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ABSTRACT 

We have gathered a sample of 112 main-sequence stars with known debris 
disks. We collected published information and performed adaptive optics obser- 
vations at Lick Observatory to determine if these debris disks are associated with 
binary or multiple stars. We discovered a previously unknown M-star companion 
to HD 1051 at a projected separation of 628 AU. We found that 25 ± 4% of our 
debris disk systems are binary or triple star systems, substantially less than the 
expected ~ 50%. The period distribution for these suggests a relative lack of sys- 
tems with 1-100 AU separations. Only a few systems have blackbody disk radii 
comparable to the binary/triple separation. Together, these two characteristics 
suggest that binaries with intermediate separations of 1-100 AU readily clear 
out their disks. We find that the fractional disk luminosity, as a proxy for disk 
mass, is generally lower for multiple systems than for single stars at any given 
i age. Hence, for a binary to possess a disk (or form planets) it must either be a 

very widely separated binary with disk particles orbiting a single star or it must 

>- be a small separation binary with a circumbinary disk. 

OO 

Subject headings: binaries: general — infrared: stars — planetary systems: for- 
mation 



j> 1. Introduction 

X 

Planet formation occurs around young stars in disks that are rich in gas and dust, 
some of which can be used to form Jovian-class planets. This needs to occur fairly rapidly, 



since disk gas generally dissipates over a period of a few million years (Haisch et al. 2001 



Uzpen et al. 2009). Eventually, even the dust in the system will be removed, either through 



accretion onto larger objects, stellar winds, or radiative processes. However, if the system 
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has formed planetesimals or larger-sized objects, collisions can occur and produce a second 
generation of dust. These dusty systems, known as debris disks, could then contain detectable 
quantities of dust with little or no gas present and they would be older than their gas-rich 



counterparts (Zuckerman 2001 Wyatt 2008, and references therein). The first debris disks 
were found by the Infrared Astronomy Satellite (IRAS). IRAS surveyed almost the entire sky 
at 12, 25, 60, and 100/iin and discovered infrared excesses around many stars including Vega, 
Fomalhaut, Pictoris, and e Eridani. The solar system's own Kuiper Belt may be analogous 
to these circumstellar disks (e.g., Luu & Jewitt 2002, and references therein). A complete 



understanding of the formation and evolution of planetary systems requires knowledge of 
the properties of disks, from the gas-rich protoplanetary disks to the gas-poor debris disks. 

Over 500 extrasolar planets have been discovered so far. Most extrasolar planet searches 
have utilized the precision radial velocity technique; however, the Kepler satellite has recently 



reported over 1200 candidate planets (see Borucki et al. 2011). While highly successful, 
the radial velocity technique is generally applied to single stars or very widely separated 
binaries. The problem is that the spectra of close binary stars are highly variable due to 
their orbital motion which typically leads to relatively large velocity uncertainties. However, 



the TATOOINE radial velocity survey (Konacki et al. 2009) is searching for planets around 



double-lined spectroscopic binaries and has reached precisions of a few m/s, comparable to, 
albeit not as good as, precisions reached for single stars. Although Konacki et al.l (2009) 



have yet to report a circumbinary planet and their sample size is small (10 systems), their 
work demonstrates that future planet searches may be performed on close binaries. 



Duquennoy & Mayor (1991) estimate that 57% of G stars may be in multiple systems. 



Other, more recent surveys show single stars are somewhat more common, but the fraction 
of multiples is still about ~ 50% (lEggleton & Tokovinin||2008| iRaghavan et al.||2010[). A 



question naturally arises: if so many stars are in multiple systems, what does this say about 
the formation and evolution of disks and planets? About 20% of known extra solar planets 



reside in wide separation binaries (Raghavan et al. 2006; Eggenberger et al. 2007), most of 



which have separations of 100s of AU. Eclipse timing variations of HW Virginis, CM Draconis, 
and NN Serpentis (a post-common envelope binary) suggest that planets may orbit these 
binaries ( Deeg et al~||2008 Lee et al]|2009 Beuermann et al.||2010 ). Despite these efforts (see 



also Konacki et al. 2009), detection of planets in circumbinary orbits remains a challenge. 



While this manuscript was being prepared, Doyle et al. (2011) presented the discovery of 



Kepler- 16b, the first circumbinary planet among the Kepler data. Both stars and the planet 
share a common orbital plane, suggesting the planet formed in a circumbinary disk. Indirect 
evidence implies the existence of rocky planets orbiting the close (3.4-day) main sequence 
binary BD +20 307. This system displays a large quantity of warm dust in the terrestrial 



planet zone ( Song et al. 2005 ) that likely is the aftermath of a collision of two rocky planets 
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that orbit this ~1 Gyr old binary system (Zuckerman et al. 2008b Weinberger et al. 2011). 



The study of circumstellar and circumbinary disks, then, can be used to comment on the 
process of planet formation around binary stars. 

There has been some previous effort to address the issue of dusty disks in binary sys- 
tems. Sub-millimeter studies of young (< 5 Myr) binaries have shown that binary stars 
with intermediate separations (1 < a < 50 — 100 AU) have lower sub-millimeter fluxes than 
more widely separated binaries or single stars (Jensen et al. 1996). Interferometric obser- 
vations of nearby ~8 Myr-old disks in a triple and quadruple system have shown the disks 



to be truncated by the nearby stellar companions (Andrews et al. 2010), as expected from 
numerical simulations (Lubow & Artymowicz 2000; Artymowicz & Lubow 1994). Among 



pre-main sequence stars, several studies have found the disk lifetimes of small to moderate 
separation binaries to be shorter than that of single stars or very widely separated binaries 



(Bouwman et al. 2006 Cieza et al. 2009). Recently, Kraus et al. (2011) presented results in 



the Taurus-Auriga star forming region demonstrating that binaries with separations <40 AU 



readily disperse protoplanetary disks. Trilling et al. (2007) used the MIPS camera on the 



Spitzer Space Telescope to search for infrared excess among 69 known binaries. They not 
only found that some binary systems have debris disks, but also that that the incidence of 
debris disks among binaries is marginally higher than for single AFGK stars older than 600 
Myr. 



2. Sample 



We have approached the question of stellar multiplicity among debris disk systems from 



a different direction than did Trilling et al. (2007). Whereas they began with a asample 



of known binarity, we instead selected a sample of stars with known infrared excesses that 
satisfy two additional criteria criteria: ages older than 10 million years, to reduce the chances 
that we include protoplanetary disks, and distances within 100 parsecs, to ensure we have 



sufficient information. We constructed our sample from Rhee et al. (2007), Rebull et al. 



(2008), Chen et al. (2005), and included BD+20 307 from Song et al. (2005). We note that 



Mawet et al. (2011) recently demonstrated that e Cephei (HIP 109857) does not contain a 
debris disk and is thus not included in our sample. 

Many Spitzer studies tend to be biased against binaries, so we were careful to select 



samples that were not biased in favor or against binaries; the Rebull et al. (2008) study was 



based on stars in the (3 Pictoris moving group, while Chen et al. (2005) searched for debris 



disks among nearby, young (12-600 Myr) stars. For the 60/zm IRAS sample from Rhee et al. 



(2007), we did not include objects marked as new candidate debris disks (see Note 2 in their 



-4- 



Table 2). By drawing from papers that use the whole sky with IRAS or in moving groups 
with Spitzer, we expect to avoid any significant bias in multiplicity fraction. 

From these references we obtain stellar properties such as spectral type, age, fractional 
infrared luminosity (Lm/L*), dust temperature, and dust orbital semi-major axis. The dust 
properties were derived in the respective references from blackbody fits to the excess emission 
after modeling the stellar photosphere (however, see Section [5]). We note that the infrared 
excesses for the sample in Trilling et al. (2007) were often faint, unlike those in our sample. 



Our efforts resulted in a catalog of 112 systems with spectral types essentially in the 
range B8 to K2, though most are A and F-type stars (see Figure [T]). Our sample, and the 
gathered information, is listed in Table [TJ The majority of these stars are IRAS detections 



from Rhee et al. (2007). Many, though, have been confirmed by Spitzer. Two of our objects 



overlap with the Trilling et al. (2007) sample: HIP 15197 and HIP 66704. Table[T]stars have 
distances ranging from 3 (e Eri) to 100 pc with a median distance of 38 pc. 



3. Procedure 



After specifying the sample, we searched the literature to determine which stars are 
known to be multiples. We used a variety of catalogs to search for information, including 



the Ninth Catalogue of Spectroscopic Binary Orbits (Pourbaix et al. 2004), the Hipparcos 



and Tycho Double Star catalog (ESA 1997), the Sixth Catalog of Orbits of Visual Binary 



Star^J the Washington Double Star catalog ( Mason et al.||2001 ), the Catalog of Components 
of Double & Multiple stars ( Dommanget fc Nys||2002 ) , the Multiple Star Catalog (Tokovinin 
1997), as well as checked SIMBAD and VizieR for papers on individual systems. While 



several systems have widely-separated candidate companions, most can be readily confirmed 
or ruled out by examining proper motions and estimated distances. Including our work in 
Section [4], we found 28 binary or triple star systems within the sample and list these multiples 
in Table |2} This corresponds to a multiplicity of 25 ± 4% (independent of whether they are 
binaries or triples), where the 4% error, as well as the uncertainties given in Table [3| are 



estimated as described by Burgasser et al. (2003). We have broken down our multiplicity 



fraction by spectral type in Table 3, which also includes the multiple fractions obtained by 
Eggleton & Tokovinin (2008) and |Duquennoy fc Mayor (1991). Our triples are hierarchical, 
with a close pair orbited by a more distant star. 

We collected period and/or semi-major axis information for our multiples. Table [2] lists 
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our separations for the binary and multiple stars. For those without periods, we estimate 
the period assuming the projected separation is the semi-major axis and the orbits are 
circular. The mass of the primary is estimated from its spectral type and age. We plot all 
the available, or derived, periods as a histogram in Figure [2] The dashed line represents the 



period distribution of Duquennoy & Mayor (1991) normalized to contain about 56 systems 



the expected number of multiples if the multiplicity fraction were 50%. This dashed line 
differs from our debris disk sample, which suggests a lack of systems with periods of about 
10 2 — 10 6 days or semi-major axes about 1-100 AU for sun-like stars. 



As mentioned in Section [Ll [Jensen et al. (1996) found that young (< 5 Myr) binary stars 



with separations 1 < a < 50 — 100 AU have lower sub-millimeter fluxes than more widely 
separated binary stars or single stars; this is approximately the same separation range as the 
gap in our sample and suggests that at such separations binary companions are effective at 
disrupting the formation of disks or accelerate the clearing of dust. However, it could also 
be that our sample is missing 20-30 binaries or multiples with intermediate-size semi-major 
axes. While most of our sample stars are well documented in the literature, there are some 
with very little information of the sort that can be used to determine their binary nature. 
We describe in the following section a search for missing companions. 



4. Adaptive Optics Search for Companions 

Wide separation companions would already have been detected with data from the 
Hipparcos satellite or proper motion surveys. To detect closer companions, one can use high- 
resolution spectroscopy over several epochs to see if the radial motion of the star changes 
indicating the presence of a massive companion. However, as most of our stars with little 
information are A-type stars, precise radial velocities are difficult to measure due to a relative 
lack of lines and rotationally-broadened profiles. An alternative way to search for companions 
is to use adaptive optics (AO). IRCAL, the infrared camera for adaptive optics at Lick 
observatory, has field of view of 20" and thus can detect companions out to nearly 1000 AU 
at 100 pc. With good seeing, companions as close as ~0'.'15 can be imaged. Because the 



peak of the binary distribution is expected to be between 10-100 AU (Duquennoy & Mayor 



1991 Eggleton & Tokovinin 2008) and our systems are located within 100 pc, IRCAL will 
be capable of finding these missing companions, should they exist. Roughly equal mass 
companions can be very quickly ruled out or confirmed after exposing for only seconds in 
the raw data alone. Fainter companions can be detected by co-adding dithered exposures. 

We had three AO runs at Lick with IRCAL (June, October 2009; August 2010). 
Additional runs were scheduled in April 2010 and March 2011, but these suffered from 



-6- 



bad weather (rain/snow). These runs were primarily intended for observation of Herschel 
DEBRIS targets. DEBRIS, or Disc Emission via a Bias-free Reconnaissance in the In- 



frared/Submillimeter, is a Herschel key project of which we are a part (see Matthews et al. 



2010). The DEBRIS project is observing 446 stars at 100 and 160-/im to look for far-infrared 



excesses indicative of cool dust in nearby star systems. There is some overlap between the 
DEBRIS sample and our own sample of 112 dusty stars and it has not been difficult to 
observe some additional Table [T] when the opportunity arises. Table [4] lists our IRCAL AO 
observations of 20 Table [T] debris disk systems. Observations were carried out using the Ks 
and BrG-2.16 filters. If a candidate was suspected, observations at other wavelengths (J, H, 
Fe II) were performed. BrG-2.16 and Fe II are narrow band filters centered on 2.167 and 
1.644/mi, respectively. They are similar to the Ks and H filters, but ~ 10 times narrower. 

Data reduction for these observations was carried out in the usual manner (flat fielding, 
dark subtraction, sky subtraction) using standard IRAF routines. Because the X and Y 
platescale are not the same, we used IDL routines to rescale the image while preserving flux 
to 76 mas/pixel in both X and Y though we later calibrate the plate scale using known 
binaries. In our search of these 20 systems we detected only 2 companions. One of these, 
HIP 35550 is a known triple system and was observed as part of the DEBRIS study. The 
other, HD 1051, we describe below. 



4.1. HD 1051AB System 

HD 1051 (HIP 1185) is a 600 Myr-old A7 star located 88.3 pc away. An infrared excess 
at 60 and lOOyum was detected around this star suggesting a debris disk with T dust = 40 K 



and i?dust = 173 AU (Rhee et al. 2007). This system was observed as part of our IRCAL 
AO observations on 2010 August 4 with the Ks filter for a total integration time of 20s 
(10 frames, each being 10 coadds of 0.2s exposures). A candidate companion, visible in the 
individual frames, was detected. 

To determine the projected separation for detected candidate companions, we calibrated 
our position angle (PA) offset and plate scale using calibrator stars from the the Sixth Catalog 
of Orbits of Visual Binary Stars. Some of these stars were targets in the DEBRIS project, 
but we used them here for calibration. We list these stars, and the results, in Table [5j We 
measured AX and AY, the difference in X and Y between the primary and the companion, 
for these stars using tasks in IRAF's daophot package and performed a least-squares fit to: 



R. A. = AAX cos 9 - BAY sin 9 
Decl. = AAX sin 9 + BAY cos 9 
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for A, B, and 9. A, B are the X, Y plate scales, respectively, and 9 corresponds to the offset 
in position angle relative to the y-axis. The IRCAL detector is aligned such that the y-axis 
on image frames is approximately North, so 9 amounts to a rotation. Our best fit values 
are: A = 72.6 mas/pix, B = 77.5 mas/pix, and 9 = 0.52 degrees. The rms differences in the 
expected position angle and separation (p) compared to the calculated values from our fit are 
adopted as our calibration uncertainties. The uncertainties in separation and position angle 
are 52 mas and 0.9 deg, respectively. These values are accurate only for the 2010 August 
observing run. 

We measure the separation of HD 1051's companion to be 7.11 ± 0.05" at 311.19 ± 
0.18 degrees East of North. Uncertainties are standard deviations of the measurements 
performed on the individual frames and do not include the calibration uncertainties. At a 
distance of 88.3 pc, 7.11" corresponds to a projected separation of 628 AU. Circular aperture 
photometry was likewise performed on the individual frames with a 2-pixel radius (0"15) 
as this aperture maximized our signal-to-noise. We measure an apparent Ks magnitude 
difference of 5.7 ± 0.1. At such separations, we could have detected companions up to 
8 magnitudes fainter than the primary at the 5a level, corresponding to ~ O.1M ; see 
Figure [3j The primary has an apparent 2MASS Ks magnitude of 6.25 ± 0.02, which implies 
the secondary has apparent magnitude of 12.0±0.1, or an absolute magnitude of 7.3±0.1. The 



system has an estimated age of 600 Myr (Rhee et al. 2007), which implies, when comparing to 
Baraffe et al. (1998) models, that the secondary has spectral type ~M3 and mass ~ O.3M . 



In order to confirm the detected object is bound to the system, and thus a true compan- 
ion, we compare the IRCAL 2010 data with existing HST NICMOS data taken during 2007 
as part of program 11157 (PI: Joseph Rhee). The HST location is 7.22±0.09" at 311.23±0.57 
degrees. As before, uncertainties are the standard deviations of multiple measurements and 
do not include systematics: ~ 0.08" in both R.A. and Deck due to the primary's location 
behind the NICMOS coronographic spot. The IRCAL and HST data are displayed in Fig- 
ure [4] and both show the companion at approximately the same location. Using the known 
distance and proper motion of the primary, as well as the measured location of the com- 
panion, we find that the secondary has not moved relative to the primary over the 3-year 
baseline. Figure [5] illustrates this by showing the measured location of the IRCAL and HST 
data and the expected motion a stationary background object would have had relative to 
the primary over this time period. 
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5. Dust Temperatures and Fractional Luminosities 



Dust temperatures, and hence disk radii or semi-major axes, are estimated in the various 
references by fitting blackbodies to the observed infrared excess. This assumes the grains 
are not much smaller than the wavelengths at which they principally emit. For objects 



with a detected excess only at 60/xm, Rhee et al. (2007) assigned a temperature of 85 K 
corresponding to the peak of a blackbody at that wavelength in F v . There are about ~20 
such objects in our sample. For objects with MlPS-measured excess emission only at 70/iin, 



Rebull et al. (2008) set a temperature of 41 K corresponding to the peak for \F X while 



Chen et al. (2005) use a temperature of 40 K, based on a modified blackbody fit to the 



dust around AU Mic. As most of our data come from Rhee et al. (2007), we have re-fit 



the spectral energy distribution (SED) for those systems with only 70/im excesses (J. Rhee, 
2011, private communication). The modified objects are HIP 2072, HIP 25486, HIP 66704, 



and HIP 92680. In comparison, Trilling et al. (2007) determine upper limits on temperature 



by using the MIPS 3a upper limit on the 24/xm emission for those objects with only 70/im 



excesses. 



For blackbody-like grains, we estimate disk semi-major axes from the stellar radius and 
temperature and the dust temperature: 



R 



dust 



R* 
T 



F 



T, 



dust 



The distributions of dust temperature and fractional luminosity are shown in Figure [6] 
While the dust temperatures are similar for both single and multiple systems, the fractional 
luminosities for single stars are, on average, 2-3 times larger than those for multiple stars. 
We applied the Kolmogorov-Smirnov (KS) test for our sample of dust temperatures and 
fractional luminosities. For the dust temperature, we find P = 0.04, which, while small, 
may not allow us to rule out that they are drawn from the same distribution. However, 
for fractional luminosity we find P = 2 x 10 -4 suggesting that the distribution of fractional 
luminosity between single and multiple stars is different. 

The difference in fractional luminosity may imply that the dust is being cleared more 
readily in multiple systems. However, it may also be that there is a difference in the age 
of single and multiple systems in our sample. Older systems are known to have disks with 



Rhee et al. 


2007 


Wyatt 


2008) 



distribution of ages in Figure [7} At a glance, the figure suggests our multiples are older, 



which is consistent with the results in Trilling et al. (2007) in that the incidence of debris 



disks around multiples is marginally higher than that for single stars at ages >600 Myr. 
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However, closer examination (the bottom panel of Figure |7j) reveals that at any given age, 
single stars have higher fractional luminosities than multiple stars. Hence, age alone cannot 
account for the lower fractional luminosities observed in multiple systems. If the incidence 
of dust around multiple stars relative to its incidence among single stars is in fact more 
likely for older stars, as suggested by our results and those of Trilling et al. (2007), then 
an explanation may reside in the long-term orbital stability of rocky objects. Studies have 
demonstrated that planetesimal orbits can be disrupted by the gravitational influence of 
planets in the system as they migrate (see, for example, Gomes et al. 2005 and references 
therein). However, as in our solar system, planetary systems should relax as they age and 
become progressively more dynamically quiescent. In contrast, the destabilizing influence of 
the gravity of stellar companions never really goes away. The very dusty system BD+20 307, 
where rocky planets orbits were likely altered following a Gyr of evolution, probably contains 
three stellar members (Zuckerman et al. 2008b). 



6. Separations 



While the period distribution is suggestive, it is instructive to compare the stellar sep- 
aration with that of the dust's orbital semi-major axis. As previously mentioned, all the 
source references for our sample calculate dust semi-major axes from fits to the spectral 
energy distribution (SED) assuming blackbody dust grains. Figure M plots the stellar sepa- 



rations and dust semi-major axes for Table 2 stars and those from Trilling et al. (2007). Only 



the shortest separation for multiples is plotted with the exception of HIP 81641, where the 
dust is located around the single primary. The grey region denotes dust locations that would 



be unstable based on the binary separation. This is adopted from Trilling et al. (2007) and is 



based on the stability study performed by Holman & Wiegert (1999). The stability of a test 



particle is defined by a critical semi-major axis, which is the maximum or minimum distance 
(depending on whether it's a circumstellar or circumbinary orbit) at which a test particle 



survives for 10 times the binary period (Holman & Wiegert 1999). The stability depends 



on the binary mass ratio and eccentricity, but we adopt average values of 0.5 for both and 



use Equations 2 and 5 in Holman & Wiegert (1999) to find that the critical semi-major axis 



is 0.12 and 3.8 times that of the binary separation. This is comparable to the values of 



0.15 and 3.5 used in Trilling et al. (2007). We note that for circular orbits, the range in 



critical semi-major axis is narrower: 0.27 and 2.3 times that of the binary separation, and as 
such we adopt the more conservative estimate with e = 0.5 in Figure [8} In this regime, the 
companion would quickly disrupt the orbits of particles be they dust or larger objects. We 
note that the actual boundaries of this unstable region are fuzzy for several reasons: (1) as 
already noted they depend on individual stellar parameters, (2) only gravitational forces are 
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considered (ie, no radiative effects), (3) mean-motion resonances can disrupt test particles 
within an ostensibly stable region, and (4) several of our systems are known triples, whereas 



the Holman & Wiegert (1999) study applies for binaries. 



Among our sample, there are 5 binaries and 1 triple that lie in the unstable region. 



In addition, 3 systems from Trilling et al. (2007) lie in that region (HD 46273, HD 80671 



HD 127726), as described in their paper. Those three systems (two are triples, one is a 
quintuple system) have detected excesses only at 70/mi (in one case only marginal). The 
ratio of dust to star separation for these 9 systems is between about 0.3 and 2.6. These 
are listed in Table |6j where the separation listed for these systems is the conflicting one. 



As Trilling et al. (2007) suggest, these systems may be undergoing a transient event, where 
dust has being generated by larger objects located farther out and has now migrated inward 
via Poynting-Robertson drag to its observed location. In addition, both the radiation and 
gravity fields for triple or higher order multiples may be sufficiently complex to affect the 
dust temperature or orbital configuration of asteroids or planetesimals. However, we describe 
below a number of effects that can potentially make systems that actually are stable instead 
appear to be in the unstable region in Figure [8] when viewed by a distant observer. 

We note that with the exception of HD 1051, the systems listed in Table [6] have infrared 
excesses detected at only 60 or 70/zm with detections at shorter wavelengths, such as 24 
and 25/im, being consistent with the stellar photosphere. Cooler dust temperatures and 
thus larger semi-major axis are consistent with the data. For clarity in Figure [8j only 
systems with a single detected excess at 70/im have right pointing arrows. Ill-determined 
dust temperature and locations are likely to be the single most important reason for systems 
to lie in the unstable zone. That is, the systems may actually not be unstable: the dust may 
be farther away and cooler than anticipated. This would shift systems towards the right 
in Figure [8} Additional observations at longer wavelengths, such as with the Stratospheric 
Observatory for Infrared Astronomy (SOFIA) and the Herschel Space Observatory, will be 
key in determining temperatures for these systems. 

An additional potential source of error in the horizontal placement of systems in Fig- 
ure [8] is the assumption of blackbody grains. Small grains radiate inefficiently and, at a 



given distance from a star, will have higher temperatures than blackbody grains (Zuckerman 



2001 and references therein). As such, at a given temperature small grains will be located 
further to the right on the Figure [8] horizontal axis than will blackbody grains. A handful of 
systems have been resolved in thermal emission (where the dust particles themselves radiate 
at wavelengths between 20 and 1000 jum), or in scattered light, where small grains reflect and 
scatter light from the central star. We consider these systems in Tableland Figure |9j noting 
that a handful of known resolved systems are not listed. These include (3 Pic, r Ceti, 7 Oph, 
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and ( 2 Ret, as either too little information is available or, in (3 Pic's case, too much, and 
thus meaningful determination of a characteristic resolved radius is not possible. Note that 
HD 141569A and HD 98800B appear in Table [7j these usually have been classified as transi- 
tion disks. For Table [7] and Figure [9] we use an average for the range in location of the dust 
or adopt the location of peak emission. Because grain semi-major axes plotted in Figure [8] 
are based on thermal emission, the relevant comparison of semi-major axes in Table [7] and 
Figure [9] is with spatially resolved thermal emission and not with scattered light. Thermally 
resolved disk radii range from one to five times larger than expected from blackbody-fit radii 
and likely are due to the particular composition, size, and porosity of grains in the disk. For 
actual disk semi-major axes 5 times larger than the blackbody model, one can shift at most 
5 systems out of 9 from the unstable region in Figure [8] Hence consideration of the actual 
sizes of disks would result in shifting a few systems to the right (larger disk radii) and into 
stability. Similar considerations could shift a few apparently stable systems into an unstable 
configuration. 

Because of projection effects and location along the orbit, the observed binary separation 
is a lower limit and could be somewhat larger than that plotted in Figure [8] The ratio between 
projected separation (p) and semi-major axis (a) for a circular orbit is: 



p = a a/ 1 — sin 2 to sin 2 i 



where to is the location along the orbit and thus has a total range of to 2ir (from Macintosh 



1994 however, see also |Leinert et aL~| [l993). The inclination, i, has a range of to tc/2 (i = 



denotes a face-on orbit). We can estimate the average ratio: 

c7r/2 p2n 
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So the average binary separation could be ~ 20% larger than what is actually measured. 
However, on an individual-case basis, objects may have inclinations and orbit locations such 
that objects move outside the unstable region. We generated 100,000 random inclinations 
between and 90° and orbit locations between and 360° to compute the ratio of projected 
to actual separation. This distribution is integrated and normalized in order to estimate how 
often ratios lower than a particular value appear. For the unstable systems listed in Table [6j 
we estimate what semi-major axis is required to yield a dust radius-to-star separation ratio of 
0.12 and calculate the probability (the final column) of achieving stability in Figure [8] based 
on the measured stellar separation. In all cases the likelihood is low that the semi-major 
axis of the binary system is much larger than the observed binary separation. Thus, at most 
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one unstable disk system would be shifted upward in Figure [8] into a stable configuration. 
However, as noted in the first paragraph of the present Section, Figure [8] assumes a binary 
eccentricity of 0.5. If the HD 1051 binary were instead to have a circular orbit, then the 
probability that the orbiting dust is stable would be 60%. In addition, should the HIP 76127 
binary also possess a circular orbit, then the circumbinary dust would be stable as its dust 
radius-to-star separation ratio is larger than the critical semi-major axis for circular orbits 
(2.3). In practice, as most binary orbits are eccentric, the companion will spend more time 
away from the primary, resulting in a lower ratio between projected separation and semi- 
major axis than the average estimated here. 

Summarizing the above considerations, data points in Figure [8] may be shifted as a 
result of incorrect temperature determination, dust grain characteristics, and projection 
effects. The dominant effect in our sample is likely the incorrect determination of dust 
temperatures, and thus, of dust semi-major axes. The other two effects are, in general, 
probably insufficient to shift all systems outside the region where dust particle orbits are 
expected to be unstable. The result of binary interactions in the disks will be to either 



to clear gaps or to truncate the disks (Lubow & Artymowicz 2000 Artymowicz & Lubow 



1994). The small number of these disk systems in the unstable zone suggests that binaries, 
although they can form disks, are more likely to disrupt the formation and evolution of 
planetary systems when their separations are comparable to typical disk sizes. 



Disk Masses 



A final comparison one can make for our debris disk sample are dust masses determined 



from (sub) millimeter photometry. Nilsson et al. (2010) have summarized all submm and 



longer wavelength detections of debris disks, but only briefly comment on binarity in their 
paper. We have compared our sample to their list, but only 23 of our debris disks have 
submm or mm measurements. Of those 23 only 3 are binaries and thus our comparisons are 
limited by small number statistics. Disk dust masses can be readily estimated with: 

S y D 2 
Md = k v B v {T d ) 

where S v is the submm flux, D is the distance between the dust and Earth, and B u is the 
blackbody flux at the dust temperature T d . The dust opacity, k u , is extrapolated from an 
opacity of 1.7 cm 2 g -1 at 880 jum assuming a dust opacity slope j3=l: k u = 1.7(880/A) /3 (see 



Zuckerman et al. 2008a). 



We list our estimated disk dust masses in Table |8j The average dust mass for our 20 
single stars is ~ 7 Moon masses (M Moon ; the median is ~3 M Mo0 n)- In contrast, our 3 
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binary star systems have an average disk mass of ~ 0.7 Mm ooi1 . These single and binary star 
dust masses, while few, are consistent with our results for fractional luminosities and dust 
separations. That is, binaries or multiples are more effective at clearing their disks and as 



such, they possess lower disk masses (as previously found for very young systems; see Jensen 



et al. 1996). We note, however, that these few binaries have ages > 200 Myr, whereas our 
single stars are predominantly younger than ~ 200 Myr. Hence, this disk mass difference, 
while suggestive, must be taken with care. 

Further submm and mm observations of our debris disk sample, for example with JCMT, 
APEX, or ALMA, will be key in gathering a larger sample of disk mass measurements for both 
single and multiple debris disk systems. In particular, for the systems in the unstable region 
of Figure |8j we extrapolate what flux density is expected at 850/im. We use the systems 



measured by Williams & Andrews (2006) and compare the -F 8 5o/-F 6 o ratio as a function of 
dust temperature in order to estimate F 85Q . We find that our unstable disks listed in Table [6] 



(not including the three Trilling et al.||2007 systems), are anticipated to have flux densities of 
~ 1 — 10 mJy at 850/im. In addition, these systems should have angular radii (from the SED 
fits) of ~ 0.4 - A". In a single hour, at 345 GHz with a 1" beam, full ALMA (54 antennae) 
can reach sensitivities of ~0.02 mJy/beam, making systems like these readily detectable and 
resolvable. 



8. Conclusions 

We find that the fraction of stars in binary or multiple systems among our debris disk 
systems is 25 ±4%. This is less than the anticipated ~ 1/2 of multiple stars and could be due 
either to a physical difference (ie, less multiples among debris disk systems) or to incomplete 
multiplicity data for some of our stars. We performed an adaptive optics search on 20 not 
previously well-studied systems in order to search for binaries in the 10-1000 AU separation 
range. We discovered only a single binary: HD 1051, an A-star with an M-star companion 
with projected separation of 628 AU. The distribution of systems in stellar separation vs. 
blackbody dust semi-major axis confirms the theoretical expectation that dust will not form 
or survive at separations comparable to that of the semi-major axis of the secondary star. 
Additionally, the fractional luminosities and disk masses, the latter inferred from a handful 
of submm detections, are lower for the multiples systems; this serves to strengthen the idea 
that these systems clear out their disks more effectively than do single stars. 

What does this imply for planet formation? First, many binaries do possess disks and as 
such could very well form planets. However, our results suggest that binaries are less likely to 
possess long-lived, massive disks. We infer this by our result of fewer than expected binaries 
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in our sample and the lower disk masses and fractional luminosities among the binaries. 



The fact that the binary semi-major axis distribution peaks around ~30 AU (Duquennoy 



& Mayor 1991) suggests that most binaries will only rarely form planets at semi-major axes 
comparable to the gas giants in our own solar system. For a binary to possess planets it 
must either be a very widely separated binary with planets orbiting a single star (in fact, 



^20% of known extrasolar planets are in such systems; see Raghavan et al.|[2006 ) or it must 



be a very closely separated binary with planets forming in a massive circumbinary disk. 

Results from unbiased Herschel surveys followed by ALMA imaging should address 
uncertainties associated with Figures [8] and [9] of the present study. 
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Table 3. Multiplicity Fractions 



Spectral Type Number Percent From ET08 From DM91 



B 


5/6 








A 


12/47 


26+^ 


46.0 




F 


5/41 




47.4 




G 


4/12 


qq+15 
°°-10 


45.0 


57 


K 


2/5 


401S 


29.1 





Note. - Fraction of multiple stars (all as percentages) broken 
down by spectral type of the primary star with comparisons from the 



literature (ET08: Eggleton fc Tokovinin 2008 DM91: Duquennoy k 



Mayor 1991). 



Table 4: IRCAL Debris Disk Observations 



Name 


UT Date 


Filter 


Exposure Time (s) 


HIP 


63076 


2009 


June 7 


Ks 


11.4 


HIP 


70952 


2009 


June 7 


J,H,Ks 


20, 11.4, 11.4 


HIP 


71284 


2009 


June 7 


BrG-2.16 


20 


HIP 


76635 


2009 


June 7 


Ks 


20 


HIP 


78554 


2009 


June 7 


Ks 


11.4 


HIP 


81800 


2009 


June 7 


Ks 


11.4 


HIP 


83480 


2009 


June 7 


Ks 


20 


HIP 


85157 


2009 


June 7 


Ks 


11.4 


HIP 


85537 


2009 


June 7 


Ks 


11.4 


HIP 


87108 


2009 


June 7 


BrG-2.16 


40 


HIP 


106741 


2009 


June 7 


Ks 


40 


HIP 


5626 


2009 


October 29 


Ks 


52 


HIP 


15197 


2009 


October 29 


BrG-2.16 


52 


HIP 


10670 


2009 


October 29 


BrG-2.16 


26 


HIP 


32480 


2009 


October 29 


BrG-2.16 


56 


HIP 


35550 


2009 


October 29 


BrG-2.16 


54 


HIP 


41152 


2009 


October 29 


Ks 


26 


HIP 


41307 


2009 


October 29 


BrG-2.16 


52 


HIP 


69732 


2010 


August 3 


J,Fe II,BrG-2.16 


6, 4.2, 3 


HIP 


1185 


2010 


August 4 


Ks 


20 
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Table 5. IRCAL Calibrations 



Name 


Grade 


PA 


P 


AX 


AY 


Calc. PA 


Calc. p 






(deg) 


(") 


(pixels) 


(pixels) 


(deg) 


(") 


HD 212698 


4 


41.4 


1.30 


-11.76 


11.98 


43.1 


1.26 


HD 165341 


1 


130.8 


5.77 


-60.70 


-47.76 


130.6 


5.76 


GJ 65 


3 


43.1 


1.96 


-18.72 


19.98 


41.8 


2.06 


HD 38 


4 


186.6 


5.94 


7.73 


-77.12 


185.9 


6.00 


HD 38 


4 


186.2 


6.04 


7.73 


-77.12 


185.9 


6.00 


HD 133640 


2 


60.1 


1.54 


-18.40 


9.89 


60.7 


1.54 


HD 160269 


3 


315.7 


0.99 


9.32 


8.73 


315.5 


0.96 


HD 146361 


4 


237.5 


7.12 


81.90 


-49.00 


237.9 


7.06 



Note. - Stars used as part of our IRCAL calibration, as described in Sec- 
tion 4J. Grade denotes the quality of the orbit, with lower numbers representing 
better orbits. The two columns after Grade denote the position angle (PA) and 
separation (p) for the system at the time of our observations. AX and AY are 
the pixel separation in the X and Y direction and the two columns following 
that are the calculated position angle and separation using the best-fit solution in 
Section El 



Table 6: Unstable debris disk systems 



Name 


Sep. 


-Rdust 


Ratio 


Type 


Stability 




(AU) 


(AU) 






(%) * 


HD 1051 


628 


173 


0.28 


binary 


6 


HIP 35550 


126 


71 


0.57 


triple 


1 


HIP 42430 


34 


21 


0.62 


binary 


1 


HIP 76127 


66 


171 


2.58 


binary 


0.1 


HIP 90185 


106 


155 


1.46 


binary 


0.2 


HIP 92680 


20 


11 


0.57 


binary 


1 


HD 46273 


26 


16 


0.62 


quintuple 


1 


HD 80671 


3.4 


2.9 


0.87 


triple 


0.6 


HD 127726 


14 


28 


1.96 


triple 


0.1 



List of unstable systems with the ratio of the dust to star separation. For completeness, we include 
the 3 systems from Trilling et al. (2007) at the bottom of the list. The final column denotes the 
probability that the inclination and location of the secondary star along a circular orbit are such 
that in actuality the system lies outside of the unstable zone (see Section [6] for details). 



-28- 



Table 7. Resolved Debris Disks 



Name 


Rrcs 


A 


Type 


Tdust 


^BB 


Refs 




(AU) 


(jaa) 




(K) 


(AU) 





HD 139664 


83 


0.6 


scattered 


80, 50 


21, 61 






Kalas et al.|2006 


Chen et 


al. 


|2009 


HD 15115 


430 


0.61 


scattered 


65 


35 






Calas et al.| 


2007a 


|Rhee et 


al. 


2007 


HD 61005 


210 


1.1 


scattered 


58 


17 


Hincs 


et al.|2007[ 


Hillenbrand et 


al. 


2008 


HD 53143 


82 


0.6 


scattered 


80 


9 






Kalas ct al. 


2006 


|Rhee et 


al. 


2007 


AU Mic 


130 


0.6 


scattered 


50 


9 






Kalas et al. 


2004 


|Rhee et 


al. 


|2007 


HD 207129 


163 


0.6 


scattered 


55 


27 






Krist et al. 


|20T0| 


|Rhee et 


al. 


2007 


HD 181327 


86 


1.1 


scattered 


75 


25 


Scl 


ncidcr et al. 


2006 


|Rhee et 


al. 


2007 


HD 15745 


300 


0.59 


scattered 


85 


22 


1 


<alas et al.l 


!Q07b| 


|Rhee et 


al. 


2007 


HR 4796A 


70 


1.1 


scattered 


110 


30 


S 


c 


neider et al. 


1999 


|Rhee et 


al. 


2007 


Fomalhaut 


141 


0.6 


scattered 


65 


73 


|Kalas et al. 


2005 


|Rhee et 


al. 


2007 


HD 141569A 


185 


1.1 


scattered 


110 


24 


Weinberger et al. 


1999 


|Rhee et 


al. 


2007 


HD 92945 


54 


0.6 


scattered 


45 


23 


Golimowski et al. 


|20TT| 


Rhce ct 


al. 


2007 


Vega 


85 


70, 160 


thermal 


80 


93 




Sibthorpe et al. 


poio 


Rhcc ct 


al. 


2007 


HD 107146 


97 


1300 


thermal 


55 


29 




ICorder et al. 


2009 


|Rhee et 


al. 


2007 


HD 32297 


115 


12-19 


thermal 


85 


28 


Mocrchen et al. 


|2007j 


|Rhee et 


al. 


2007 


HD 10647 


85 


70 


thermal 


65 


22 






Liseau et al. 


2010 


Rhee et 


al. 


2007 


HR 8799 


200 


70 


thermal 


150, 45 


9, 95 












ISuet 


al. 


2009 


HD 191089 


59 


18.2 


thermal 


95 


15 


Churcher et al. 


|201l||Rnee et 


al. 


2007 


jj Tel 


24 


18.3 


thermal 


150 


15 






Smith et al. 


|2009| [Rheeef 


al. 


2007 


e Eri 


60 


850 


thermal 


40 


27 




G 


reaves et al. 


lggsHRheeef 
1 i W— 


al. 


2007 


r\ Crv 


145 


100 


thermal 


354, 31 


174, 1.4 










Matthews et 


al. 


2010 


P UMa 


47 


100 


thermal 


109 


51 










Matthews et 


al. 


2010 


P Leo 


39 


100 


thermal 


112 


23 










Matthews et 


al. 


2010 


HD 98800B 


13 


880 


thermal 


160 


2.2 




Andrews et al. 


2010 


| |Low et 


al. 


|2005 



Note. — Disk radii (R res ) for resolved systems. This table is not meant to represent the most descriptive 
values for the system (for example, some disks are clearly asymmetric); we refer the reader to the individual 
references for details. For systems resolved in both thermal and scattered emission, we list only the thermally 
resolved radius. Resolved and blackbody-fit radii are compared in Figure [9] 
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Table 8: Debris disk dust masses 



Name 


\ 

'Must 


TP 


Mass 


Age 




(//m) 


(mJy) 


( A/T \ 

\iv±moon ) 


/ T\ 1 T ™\ 

(Myr) 


ATT IV /T4 „ 

AU Mic 


850 


1/1/1 _l_ 1 O 


n c 

0.5 


1 O 

12 


p Pic 


1300 


o a n i a < ■ 

24.9 ± 2.6 


r o 

5.8 


12 


Hlr 952 


o7(J 


51.7 ± 0.2 


OP cr 


1 O 

12 


TUTU) HH07Q 

Hir yyzio 


orn 

o5U 


E /I _1_ 1 IX 

54 ± 15 


2.2 


on 


TUTU) 1 1 /iion 


orn 

o5U 


1 n q _i_ i o 
1U.O ± 1.0 


P o 
O.O 


on 
oU 


TTTT) 1 p A Af\ 

Hlr lo449 


870 


17.6 ± 8 


no i 

33.1 


50 


HIP 10670 


850 


r r i "i o 

5.5 ±1.8 


1.8 


100 


TTTTD QCO/1C 


OOU 


yo ± rz 


o c; 
z.o 


i nn 
ruu 


TTTT ) 1 i o^n 

rllr llobU 


o5U 


4.9 ± l.b 


2.9 


1UU 


TTTT3 enn7/l 

Hlr oUU 1 4 


q t^n 
o5U 


on _i_ /i 
zU ± 4 


5. ( 


i nn 
1UU 


HIP 22226 


870 


6.9 ± 5 


13.4 


100 


HIP 87108 


870 


12.8 ±5.2 


2.6 


200 


HIP 101612 


870 


13 ±7.1 


3.2 


200 


HIP 6878 


1200 


3.2 ±0.9 


4.5 


200 


HIP 90936 


870 


18 ±5.4 


9.7 


200 


Vega 


1300 


11.4 ± 1.7 


0.6 


220 


T] Crv 


850 


7.5 ± 1.2 


0.3 


300 


HIP 32480 


850 


5.5 ± 1.1 


0.5 


600 


e Eri 


1300 


24.2 ±3.4 


0.4 


730 


cr Boo 


850 


6.2 ± 1.7 


0.7 


1000 


Fomalhaut 


850 


97 ±5 


1.7 


220 


HIP 107649 


870 


5.1 ±2.7 


0.5 


600 


HIP 27072 


850 


2.4 ±1 


0.04 


1000 



Disk dust masses for Table[T]stars listed in Nilsson et al. (2010 ). The last three systems are binaries, 
^dust is the wavelength at which the dust is detected and the one used to estimate dust mass (see 
Section [7]) ; dust temperatures and distances from Earth are listed Table [lj 
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Fig. 1. — Spectral type distribution of our 112-star debris disk sample. 
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Fig. 2. — Period distribution of our debris disk sample. For triple systems, we include both 
periods (for example, the A-B period and the AB-C period). For comparison, the semi-major 
axis for a sun-like star is displayed on top. The dashed line illustrates the expected period 



distribution from Duquennoy & Mayor (1991) (and very simlar to Eggleton & Tokovinin 



2008) normalized to contain the expected number of binaries or multiples in the sample 
(112 x 0.5 = 56 multiples). While the short and long period binary distribution approxi- 
mately matches the Duquennoy & Mayor (1991) distribution, there is a lack of systems with 
intermediate periods (separations ~ 1 — 100 AU). 
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Fig. 3. — Magnitude limits (5a) in the Ks filter as a function of separation for our HD 1051 
(HIP 1185) data. Limits are obtained by performing 2-pixel radius aperture photometry (the 
same as for HD 1051B) at random locations in annuli around the primary. The standard 
deviation of these measurements provides an estimate of the noise. The detected companion 
is plotted circle. 




Fig. 4.— Lick IRCAL (left) and HST NICMOS image (right) for HD 1051. 
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Fig. 5.— Measured offsets for HD 1051B in Lick IRCAL and HST NICMOS data. The black 
line illustrates the motion relative to HD 1051 of a stationary background object over the 
3-year period. The circle at the end of the black line is the predicted location of a background 
object in 2010. 
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Fig. 6. — Dust temperatures (top) and fractional luminosities (bottom) for our debris disk 
sample. Open bars represent single systems, while the shaded bars represent multiples. The 
large number of systems in the ~ 70 — 90 K bin is partially an artifact of the analysis and 
includes systems with detections at only 60 or 70/mi, see discussion in Section |5j While the 
dust temperatures are similar for both single and multiple systems, the fractional luminosities 
for single stars are, on average, 2-3 times larger than those for multiple stars. BD +20 307 
is the rightmost system in both panels (Weinberger et al. 2011). 
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Fig. 7. — Age distribution in our debris disk sample. Open bars represent single systems, 
while the shaded bars represent multiples. While the histogram (top) suggests more multiple 
systems are older compared to the single stars, the plot of fractional luminosity vs. age 
(bottom; circles are single stars, squares are multiple systems) shows that at any given age, 
multiple stars are more likely to have lower fractional luminosities compared to single stars. 
The ~1 Gyr-old binary system with Lir/L* ~ 0.03 is BD +20 307. 
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Stellar separation versus blackbody dust semi-major axis for debris disk stars 



including those from Trilling et al. (2007). Triple systems are plotted at both orbital separa- 
tions with a dotted line connecting them; the exception is HIP 81641 where the disk may be 
around the single primary so the separation of the companion close binary is not relevant. 
The two systems with right pointing arrows are HIP 66704 and HIP 92680 (PZ Tel), for 
which dust emission is detected only at 70/zm; for clarity we do not label the systems with 
excesses detected only at 60/zm (see Section [5]). The grey region in the plot highlights where 
the gravitational field of a companion is expected to significantly affect dust in the system 
(see the first paragraph of Section [6J 
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Fig. 9.— Resolved vs black body disk radii for debris disks as listed in Table [7} The 
lines represent ratios of 1:1 (solid), 2:1 (dashed), 5:1 (dotted), and 15:1 (triple-dot-dashed) 
when comparing resolved disk sizes to those inferred from blackbody fits. Circles are disks 
resolved in thermal emission, while squares are disks resolved in scattered light. Thermally 
resolved disks have sizes ~ 1 — 5 times that estimated from the blackbody SED, whereas 
disks resolved in scattered light tend to be more extended. 



